Thermodynamics and Crossover Phenomena in the Correlation Lengths of the 

One-Dimensional t-J Model 



J. SirkerQ and A. Kliimpei]^ 

Theoretische Physik I, Universitdt Dortmund, Otto-Hahn-Str. 4, Dortmund, Germany 

(Dated: February 1, 2008) 

We investigate the thermodynamics of the one-dimensional t-J model using transfer matrix renor- 
malization group (TMRG) algorithms and present results for quantities like particle number, specific 
heat, spin susceptibility and compressibility. Based on these results we confirm a phase diagram 
consisting of a Tomonaga-Luttinger liquid (TLL) phase for small J/t and a phase separated state 
for J/t large. Close to phase separation we find a spin-gap (Luther-Emery) phase at low densities 
consistent with predictions by other studies. At the supersymmetric point we compare our results 
with exact results from the Bethe ansatz and find excellent agreement. In particular we focus on 
the calculation of correlation lengths and static correlation functions and study the crossover from 
the non-universal high T lattice into the quantum critical regime. At the supersymmetric point we 
compare in detail with predictions by conformal field theory (CFT) and TLL theory and show the 
importance of logarithmic corrections. 
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I. INTRODUCTION 

The t-J model is one of the most fundamental systems 
of strongly correlated electrons. The two-dimensional 
version has attracted much attention because it is be- 
lieved that it describes the basic interactions in the 
copper-oxygen planes of high- Tc superconductors. For 
the one-dimensional (ID) t-J model much progress has 
been achieved^ |UsiHfi |-,j2arious analytical and numeri- 
cal techniques Jj'aaaiJ'ElElEl At the supersymmetric point 
J/t = 2 the model is solvable by the Bethe ansatz and 
ground state properties as weU as the excitation spec- 
tra have been obtained exactly.El Because the two critical 
excitations of spin and charge type are separated, the 
properties can be described by two independent c = 1 
Virasoro algebras. By a combination of finite-size results 
from the Bethe ansatz and conformal field theory (CFT) 
it is therefore also possible to calculate the critical ex^ 
ponents of algebraically decaying correlation functions .□ 
This explicitly shows that the t-J model at the super- 
symmetric point behaves as a Tomonaga-Luttinger liquid 
(TLL) for all electron densities. Thermodynamic quan- 
tities at this special poiat have been obtained by ther- 
modynamic Bethe ansataa as well as by a combination of 
a Trotter-Suzuki mapping leading to a, quantum transfer 
matrix (QTM) and the Bethe ansatz.tl Exact results are 
also available in the limit J/t ^ where the t-J model is 
equivalent to the Hubbatd model with U/t ^ oo, show- 
ing again TLL behavior £j It is therefore believed that the 
t-J model shows TLL properties for all < J/t <^|jsLhat 
is supported by various numerical calculations.Eialjij In 
these numerical works there is also general agreement 
that the t-J model phase separates for J/t — 2.8 to 3.5 
depending on the electron density. 

Already Ogata et al^ conjectured a third phase with 
a spin gap in the low-density region for J/t > 2. How- 
ever, by calculating the spin susceptibility for an elec- 
tron density n = 1/3 on small chains no evidence for 



a spin gap was found. Also the variational quanturp. 
Monte Carlo (QMC) calculations of Hellberg and MeleS 
could not confirm the appearance of a spin gap. By us- 
ing the same method but other trial wave functions a 
phase with Lu.ther-Emery (LE) properties was found by 
Chen and LeeQ and Kobayashi et alE at low densities and 
2 < J/t < 3.1. However, the obtained result strongly de- 
pends on the trial wave functions used in the calculations. 
Completely different phase boundaries with the spin-gap 
phase extending into the high, density region have been 
obtained by Nakamura et aln using a renormalization 
group treatment of the Tomonaga-Luttinger model un- 
der the assumption that the spin gap is caused by an 
attractive backward scattering process. They argue that 
the spin gap was underestimated in the numerical cal- 
culations, because it is the result of a marginal operator 
leading to an exponentially small gap. 

In this paper we want to examine thermodynamics 
of the one-dimensional t-J model in the whole J/t- 
parameter region. In particular, we are interested in 
the temperature dependence of correlation lengths. We 
use a grandcanonical description of the model, where the 
Hamiltonian is given by 

H = -tJ2PicU<^^:+l,.+cU,^aC.,.)P (1) 
i i i 

with a magnetic field h, a chemical potential and P 
being the projection operator onto the Hilbert-subspace 
without double occupancy. To study thermodynamic 
properties also away from the supersymmetric point the 
TMRG provides a powerful numerical tool. This method 
is particularly suited, because the thermodynamic limit 
is performed exactly and it has been appliecL-aucp^fe 
fully to various one-dimensional systems before. EdOllj'La 
Compared with QMC methods it has further the advan- 
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tage of never suffering under the fermion sign problem 
restricting QMC often to relatively high temperatures. 

In Sec. H we give a brief introduction into the calcu- 
lation of thermodynamic quantities by the TMRG algo- 
rithm and point to a novel algorithm with several ad- 
vantages. Our res ults for the supersymmetric point are 
shown in Sec. Ill where we also compare the Jiumerics 
with the Bethe ansatz results by Jiittner et ala In par- 
ticular we calculate for this special point the temperature 
dependence of various correlation lengths and static cor- 
relation functions and discuss the low-temperature be- 
havior in comparison to CFT predictions by Kawakami 
and Yangp in detail. Because the t-J model is realized 
by the Hubbard model in the limit U ^ t the case 
J ^ 2t'^/\U\ < t is physically very relevant. Results for 
J/t — 0.35, a value often used in literature, are given in 
Sec. 1^ As already mentioned the t-J model phase sep- 
arates for J/t large. The meaning of this in the grand- 
canonical ensemble is explained in Sec. |^. Without using 
any assumption about the ground state or the low energy 
effective theory, the existence of a LE phase is proven in 



Sec. VI by calculating directly spin susceptibilities as well 
as spin-spin and density-density correlation lengths. In 



Sec. VII we study the t-J model with an additional Ising- 
like anisotropy. The final section is devoted to our con- 
clusions and in appendix ^ we compare numerical and 
exact results for the density-density correlation function 
of free spinless fermions. 



II. TMRG 

After a decomposition of the Hamiltonian H into even 
(He) and odd parts (Ho) the partition function is ex- 
pressed by means of the Trotter formula 



Z = Ttc-^" = lim 



Tr 



} (2) 



with e — P/M and (3 being the inverse temperature 
leading to a classical model on a lattice with checker- 
board structure. The column-to-column transfer matrix 
(QTM),rj\/, is a non-symmetric matrix describing the 
evolution along the spatial direction. The thermody- 
namic limit with fixed Trotter number M is performed 
exactly, because the free energy of the infinite chain is 
given solely by the largest eigenvalue Aq of this QTM 



T 



■In An 



(3) 



where Aq is unique and a real, positive number for all 
temperatures. Since a vanishing gap between leading and 
next-leading eigenvalue indicates a phase transition, such 
a degeneracy is not possible for a ID quantum system at 
finite temperature. The calculations are simplified by 
the conservation laws for spin and particle number or 
equivalently for the number of particles with spin up (N^) 
and spin down {N^) leading to a block structure of Tm- 
Expectation values of local operators can be expressed in 



terms of the left and right eigenvectors belonging to the 
largest eigenvalue and calculated directly. This is used 
for the inner energy u and the magnetization m. The 
behavior of a two-point correlation function is given by 
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with Ma being matrixelements. The correlation lengths 
can be calculated from the QTM by 



and the wavevectors ka by 



Ao 

Aa 



(5) 



fca = ^ arg nTT (n = or n = 1), (6) 

where A^ are eigenvalues in the block of the QTM with 
non- vanishing matrixelements Ma ■ This condition is con- 
trolled by a direct evaluation of 



Ma = 



{^^\TM{Ao)\^^){^^\TM{Ar)\^^) 
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where {'^'q 



are the left and right eigenvectors be- 
longing to the largest eigenvalue Aq and (^'^l, 1^'^) the 
left and right eigenvectors belonging to A^ . The value for 
the wavevector ka is not unique, because in the checker- 
board decomposition of the partition function a local 
transfer matrix covers 2 sites. To overcome this ambigu- 
ity we have also applied a different Trotter-Suzuki map- 
ping leading to a classical lattice with alternating rows. 
Because the transfer matrix in this formalism can be for- 
mulated for a single column, the wavevector in Eq. (|^) 
is determined unambiguously. Details of thi£L. modified 
TMRG algorithm will be published elsewhere.Ej We have 
convinced ourselves that the accuracy of both algorithms 
is of the same order and use both of them in the follow- 
ing. During all calculations we use the infinite DMRG 
algorithm to increase the Trotter number M being equiv- 
alent to a decrease of the temperature T. The number of 
states kept in the DMRG varies between 64 and 240 and 
the fixed parameter e between 0.025 and 0.05. 



III. THE SUPERSYMMETRIC POINT 

To check the numerics we have calculated several 
thermodynamic quantities for the supersymmetric point 
J/t = 2. and compared with the exact results by Bethe 
ansatz.El We have chosen three chemical potentials corre- 
sponding in the low-temperature limit to a high, medium 
and low electron density (see Fig. |l|). Obviously the par- 
ticle density n at a given chemical potential depends on 
temperature. Counting the degrees of freedom per lat- 
tice site immediately implies n = 2/3 for any finite fj, 
in the limit T — > oo, whereas the density at finite tem- 
perature is given by n = n{p,,T). For the same chemical 
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FIG. 1: Temperature dependence of the density for three dif- 
ferent chemical potentials. The lines are given by the TMRG 
results, whereas the symbols denote the results from the Bethe 
ansatz. 
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FIG. 2: Deviation of the free energies calculated by means of 
TMRG in comparison to the exact results for different tem- 
peratures. The lines are guide to the eye. The inset shows the 
Bethe ansatz results (symbols) and the TMRG results (lines) . 



are equal to 1. According to Eq. (^) we have fitted 
the numerical data and tried to determine the errors of 
the fit parameters by a variation of the fit region from 
T/t e [0.01 : 0.02] to T/t G [0.01 : 0.05]. The estimates 
for the ground-state energies eo coincide with the exact 
results with deviations that are slightly larger than the 
errors calculated by the fit procedure (see table ||) . The 
remaining deviation is due to systematic errors. The pa- 
rameter l/vg -|- \/vc is difficult to obtain from such a 
fit, because the extent of the low-temperature region in 
which Eq. (H) is valid is very small (see Bethe ansatz re- 
sults in FigT]^) making it necessary to restrict the fit to 
the temperature region defined above. In such a small 
interval, however, a change of the parameter l/vg + 1/vc 
by a factor of 2 corresponds to deviations in the fitted 
free energy of the order 10"'^ only, so that the TMRG 
data are not accurate enough to determine this param- 
eter. However, it is possible to calculate the velocities 
from TMRG by using the results for the susceptibilities 
and correlation lengths to be discussed later on. We do 
not want to pursue this further, because the estimation 
of Vs and Vc is not our main goal. 

To calculate the specific heat c„, we first have to cal- 
culate the entropy S. This can be done by using the 
relation S = {u — f)/T, where the inner energy u is 
directly calculated in the TMRG algorithm as a local ex- 
pectation value or alternatively, by the numerical deriva- 
tive S = —df/dT. At first sight it seems better to cal- 
culate the entropy without using numerical derivatives, 
however, a local expectation value directly involves the 
eigenvectors of the QTM which are generally less accu- 
rate than the corresponding eigenvalues. By comparing 
both results with the exact ones we have convinced our- 
selves that the numerical derivative leads to more accu- 
rate results. Because in the calculations only the chem- 
ical potential fj, but not the particle density n can be 
fixed, we calculate the specific heat by using the thermo- 
dynamic relation 



potentials with = 100 states retained in the DMRG al- 
gorithm and e = 0.05 the free energies and their accuracy 
are shown in Fig. ^ in ppmparison to the Bethe ansatz 
results by Jiittner et alB In all three cases the accuracy 
is of the order 10"'' if T/t > 0.1. For lower temperatures 
0.01 < T/t < 0.1 the errors for the medium and high 
electron density grow up to 10""^, whereas the error in 
the low density case remains of the same order. More ac- 
curate results can be obtained if more states are retained 
or an extrapolation in the number of states N and the 
small parameter e is made. However, the obtained accu- 
racy is sufficient for our purposes. From CFT it is known 
that the low-temperature asymptotics is given by 

f = eo-^('- + ^)T' (8) 

6 \Vs VcJ 

where Vs^c are the velocities of the spinon and holon 
excitations, respectively. Here the central charges Cs,c 



UTj„-' [Urj, [otJ^W)^ ■ 

(9) 

The second term is again evaluated using numerical 
derivatives. This procedure enhances the numerical er- 
rors, but we are able to reproduce the exact results 
within errors of the order lO"'^ down to temperatures 
of T/t « 0.2 as shown in Fig. |. The CFT result for 
the free energy, Eq. (^) , leads to a linear temperature de- 
pendence of the specific heat at low temperatures with a 
coefficient given by tt{1/vc + l/vs)/i- The Bethe ansatz 
results clearly show this linear behavior at very low tem- 
peratures but within the TMRG this temperature region 
is not accessible because of the errors caused by the nu- 
merical derivatives used in Eq. (^). 

To calculate the spin susceptibility Xs, a small mag- 
netic field h — 10~^ is applied. From the resulting 
magnetization m the susceptibility at vanishing magnetic 
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TABLE L Velocities «s,c and ground-state energies eo from the Bethe ansatz in comparison to values from a fit of the immerical 
data. The errorbars of eo and 1/vs + l/nc correspond to the described variation of the fit region, but the errors for l/u^ + 1/wc 
are in fact much larger (see explanation in the main text). 
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FIG. 3: Specific heats as calculated by means of TMRG and 
Eq. (^) (lines) and exact results by Bethe ansatz (symbols). 
Note the two peak structure caused by spin-charge separation 
as explained in the text. 



field is evaluated by 



Xs 



h=0 



m 



Similarly the charge susceptibility Xc 
is given by 



Xc= ^ 



dn 
9/i ' 



(10) 

or compressibility 
(11) 



where again numerical derivatives are used and the vari- 
ation of ^ is typically of the order 10"^. Here we are able 
to reproduce the exact results down to temperatures of 
T/t — 0.01 as shown in Fig. ^. The high temperature 
asymptotics of the spin susceptibility is given by the pure 
paramagnetic part Xs{T oo) ~ 2s(s -f l)/9r what is 
easily understood when Xs is expressed as a sum of two- 
point spin correlations. For T = the spin susceptibility 
has a finite value as expected from the linear dispersion 
of the spinon excitations. Within CFT the zero tem- 
perature value is given explicitly by Xs = l/(27rws) in 
good agreement with the numerical results. The max- 
ima at finite temperatures are determined by the band 
structure. The charge susceptibility is given in the same 
way as sum over two-point density-density correlations 
leading to the high-temperature asymptotics Xc ~ 2/9T. 
For T = the charge susceptibility is again given by 
CFT as Xc = il{Q)l{'^'"c) with £,c{Q) being the dressed 
charge. Note, that the charge susceptibility is diverging 



H=-0.7 (exact) 
H=-1.4 (exact) 
H=-1.9 (exact) 
H=-0.7 (TMRG) 
fi=-1.4 (TMRG) 
H=-1.9 (TMRG) 




FIG. 4: The main figure shows the spin susceptibilities, the in- 
set the charge susceptibilities (compressibilities) , where again 
lines denote the numerical results and symbols the exact ones. 
The symbols at T = denote in both graphs the CFT results. 



for T ^ in the two limiting cases n ^ and n 1, 
because Dj, ~* in both cases. The spin susceptibility Xs 
shows different behavior with divergence for T ^ only 
in the limiting case n — > as only here Vs ^ Q. 

To summarize our numerical findings for c, Xs and Xc in 
a qualitative manner we may use a picture of the ground 
state as a "liquid" consisting of bound singlet pairs of 
electrons. Above this ground state there are two differ- 
ent elementary excitations. First, a momentum transfer 
onto any individual pair is possible by keeping the bound 
pair intact (holon excitation) . Second, a breaking of any 
bound pair into its components (electrons) with the two 
spin-1/2 objects coupling either to 5 = or 1 (spinon 
excitation). Once the "free" electrons are produced by 
spinon excitations, they may acquire momentum and en- 
ergy individually which we refer to as incoherent single 
particle motion. Basically, any structure observed in the 
T-dependence of c, Xs and Xc can be attributed to the 
saturation of a particular type of excitation. In Xs (Xc) 
we see a finite temperature maximum at a temperature 
of the order of Vs (vc)- Note that the characteristic holon 
temperature is always lower than or equal to the spinon 
temperature. This is rather natural as the saturation of 
spinous is marked by the practical absence of pairs such 
that saturation of holon excitations must have occurred 
at lower or equal temperature. In principle, these struc- 
tures are also visible in the specific heat. In most cases. 
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however, only the spin-maximum is clearly noticeable. 
For generic densities the small holon-maximum disap- 
pears in the low temperature regime of c{T) with rather 
steep slope due to the spinon-excitations. At lower densi- 
ties the holon and spinon structures are located at similar 
temperatures. In all cases, the higher temperature region 
is dominated by a broad maximum due to incoherent sin- 
gle particle motion. The corresponding energies are of 
the order of spinon excitations at particle densities close 
to 1 in which case the structures merge. 



A. Correlation lengths 
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spin-spin (k incomm.) 



For the correlation lengths of the supersymmetric t- 
J model no exact results are available. To obtain nev- 
ertheless a measure for the accuracy of the numerical 
data we have calculated the leading correlation length 
for a system of free spinless fermions (see appendix |A|). 
The comparison with the exact result shows that TMRG 
yields reliable results down to temperatures T/t ~ 0.1. 
Because the free fermion model is harder to tackle within 
the DMRG scheme due to the rapidly decaying spectrum 
of the reduced density matrix, we believe that for the t-J 
model a similar or even better accuracy is obtained. 

We concentrate now on the calculation of correlation 
lengths for the supersymmetric t-J model to study the 
crossover from the high T lattice into the quantum critical 
regime determined by T <C i. In the high T lattice region 
we expect non-universal properties dependent on the mi- 
croscopic Hamiltonian H, whereas the quantum critical 
regime should show universal TLL properties. Respect- 
ing the selection rules, the density-density (d-d CL's) and 
longitudinal spin-spin correlation lengths (s-s CL's) are 
in the block of the QTM with unchanged quantum num- 
bers (AA^i = ANi =0). To distinguish between them, 
the matrixelement Ma in Eq. (0) has to be calculated 
explicitly. With zero magnetic field and isotropic spin in- 
teractions, it turns out that Ma is either non-zero for the 
density or for the operator so that all d-d and s-s CL's 
are different from each other. This is a consequence of the 
SU{2) spin symmetry. If a magnetic field is applied this 
symmetry is broken and less stringent selection rules ap- 
ply such that all eigenvalues discussed above contribute 
to the s-s as well as to the d-d correlation function. On 
the other hand the singlet pair operator — C|,i+iC|_i 
changes the quantum numbers by AA^| = AA^j = ±1 and 
the eigenvalues corresponding to the singlet pair CL's are 
found in the block with these quantum numbers. The 
leading d-d and s-s CL's (times temperature) and the 
corresponding wavevectors for J/t — 2.0, /i = —1.4 and 
zero magnetic field are shown in Fig. ^. At high temper- 
atures the largest s-s CL is given by a real negative eigen- 
value leading to 7r-oscillations. However, at a well defined 
crossover temperature T,. « 0.8 a real, positive eigenvalue 
becomes largest, which is associated with fc = 0. Regard- 
ing only the largest spin CL this means that there is a 
non-analyticity at Tc. We want to point out that any 



FIG. 5: Temperature dependence of the leading d-d and s-s 
CL's for J/t — 2.0 and fi = —1.4. The triangle up (star) gives 
the zero temperature result from CFT for the non-oscillating 
s-s (d-d) and the square that for the ikp part of the d-d. The 
triangle down at zero temperature denotes the CFT value for 
the 2kF part of s-s and d-d, whereas the triangles down at 
T/t = 0.1 are given by CFT plus logarithmic corrections as 
described in the text. The inset shows the wavevectors k in 
the case of incommensurate oscillations. The circles denote 
the values for k — ±2fcF at zero temperature as expected from 
CFT. 



thermodynamic quantity derived from the free energy is 
an analytic function at finite T. Phase transitions and 
corresponding singularities only occur at T = 0. How- 
ever, quantities describing the asymptotics of correlation 
functions (i.e. CL) may show non-analyticities even at 
finite temperature. Such crossovers between CL's are 
characteristic for the non-universal high T lattice regime, 
whereas no crossovers are expected in the universal quan- 
tum critical regime described by CFT. The second largest 
spin CL is given at low temperatures by a pair of complex 
conjugated eigenvalues. Complex eigenvalues always ap- 
pear in pairs, because the QTM is nonsymmetric but real 
and lead to incommensurate oscillations, i.e. fc 0, tt. 
An incommensurate wavevector k is depending on tem- 
perature as seen in the inset of Fig. ||. The leading d-d CL 
shows no oscillations at high temperatures but is crossed 
at lower temperatures by a real, negative eigenvalue lead- 
ing to TT-oscillations before again the non-oscillating part 
dominates. Below Tc « 0.5 the largest d-d CL is given by 
a complex eigenvalue pair. The wavevector for this in- 
commensurate part is again shown in the inset of Fig. |^. 
It is illuminating to regard first the d-d correlation for 
a system of free spinless fermions. For such a sys- 
tem the CL's can be easily calculated (see appendix 
^) and it turns out that a non-oscillating part and 
a part with k = n arise, both with a CL given by 

— 2arcsinh(7rT). Therefore the CL's at high temper- 
ature behave as ~ lii(2^) consistent with the numer- 
ical results. In the low temperature limit it follows that 

^ l/(27rr) what is expected from CFT when the ve- 
locity V and the scaling dimension x are equal to 1. The 
matrixelement for the CL is given by M = 2/(7r^ + /3^) 
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indicating that M ^ when T 0. This is reasonable, 
because at zero temperature the CL's diverge and the 
matrixelements have to approach zero so that the series 
in Eq. is not diverging but sums up to an algebraic 
function. The numerical results also show that the ma- 
trixelements Ma — + if T ^ 0. 

To understand the behavior of the CL's at low temper- 
atures in detail it is useful to regard the results from the 
Bethe ansatz and theJinite-size scaling technique in CFT 
at zero temperature.Q In general, a two-point correlation 
function of the scaling fields 4>a± (^j ''') with conformal 
weights is given by 



approximately nT—,o ~ 0.524 and therefore kp « 0.823. 
In the inset of Fig. |5| the circles denote k — ±2kp at 
zero temperature indicating that the leading d-d and s-s 
CL's with incommensurate oscill atio ns correspond to the 
2fci?-oscillating part. From Eq. (12a) we derive 



(0A±(r-,T)0A±(O,O)) 



(12a) 



J{2Tr-2kF-f-2kFi)DaX„i{2TT-2kF';)DsX 



(r — Wct)'^^^ {r + Wct)'^'^': {r — Wst)^^^ (r -|- Wst)'^'^^ 

where the conformal weights are determined by quantum 
numbers of the elementary excitations. Considering here 
only the case of vanishing magnetic field and using the 
notation of Ref. ^ the conformal weights reduce to a 
simple form 



(12b) 
(12c) 



(,c{Q) is the dressed charge which can be calcu- 
lated explicitly by solving an integral equation and 
(Ic, Is, Dc, Ds, Nf., Ns) are quantum numbers of the ex- 
citation. To be specific, Ic counts the total number of 
holes. Is the number of holes with respect to the up spins, 
Dc (Ds) the number of holons (spinous) being transfered 
from one Fermi point to the other and Nc (Ns) the num- 
ber of charge (spin) particle-hole excitations. The Fermi 
momentum kp^-^ for up (down) spin electrons is given 
by 



(12d) 



where m is the magnetization. If no magnetic field is 
applied then kp-^ = kpi = kp. As shown in Fig. |l| the 
density in the zero temperature limit for ^ = —1.4 is 



(n(r)n(O)) 



const, -f Aqt ^ + cos{2kpr) 
-\A^r-<^- coa{Akpr) (13a) 



for the equal-time d-d correlation, where Ai are matrix- 
elements. The non-oscillating part is due to the low- 
est particle-hole excitation (0,0,0,0,1,0), the 2kp part 
due to a (0, 0, ±1, =f1, 0, 0) excitation, while the Akp part 
arises from the excitation (0, 0, ±1, 0, 0, 0). Thus, the 



critical exponents can be calculated from Eq. (12b, 12c) 
leading to 



1 



edQ) 



Pc = 2eAQ) 



(13b) 



The equ al-time s-s correlation has the same form as 
Eq. (13a), but the constant as well as the Akp part are 
absent and the matrixelements are different. The criti- 
cal exponent of the 2kp part is the same bui the corre- 
sponding excitation is now (0, 0, 0, ±1, 0, 0)£3 As shown 
in Ref. ^ the dressed charge ^dQ) varies between \/2 for 
n = and 1 at half-filling. Thus, ac is always the small- 
est critical exponent and therefore the 2kp parts of the 
correlation functions dominate at zero temperature. 

The results for the correlation functions shown here are 
absolutely consistent with TLL theory when the identifi- 
cation 



cm 



(14) 



is used. For the Tomonaga-Luttinger model it is also pos- 
sible to calculate multiplic ative logarithmic corrections to 
the algebraic terms of Eq. ( |l3a|) which are not directly ac- 
cessible within Bethe ansatz and CFT. It turns out that 
the 2kp parts of the d-d and the s-s correlation have dif- 
ferent logarithmic corrections, given by ln-3/2 r (In^/^ r) 
for the d-d (s-s) correlation.Ej 

The long-distance asymptotics of correlation functions 
at small finite temperature can still be obtained from con- 
formal invariance by the usual mapping of the complex 
plane onto a strip of width 1 /T with periodic boundary 
conditions. Eq. (|12al) is then replaced by 



^i{2Tr-2kF\—2kFi)DaX^i{2Tr-2kFy)D,x 



sinh(7rT(a; — iVcT)/vc) J \VcSmh{T:T{x + Wct)/vc) 



(15) 



Vs siuh{TrT(x — Wst) /vs) J \VsSmh(TrT{x + ivsT)/vs) 

I 

Thus, all CL's diverge in the low temperature limit as where the coefficient 7 — Vs/2Tr (7 — Uc/2vr) for the non- 

i = V = ^ (16) 



7 



oscillating part of the s-s (d-d) CL, 

7 = 



for both 2kp parts whereas 

Vr 

1 



(17a) 



(17b) 



for the Akp part of the d-d CL. Contrary to zero temper- 
ature the answer to the question which part dominates 
depends not only on the scaling dimensions Xs,c but also 
on the velocities Vs ,A In Fig. ||the values for these 7's are 
denoted by symbols at T = showing good agreement for 
the non-oscillating parts between the numerics and the 
CFT results. The 4fcF-oscillating d-d CL is so s mall that 
it is difficult to obtain numerically. From Eq. (17a) the 
2kF s-s and d-d CL's are expected to be equal in the low- 
temperature limit. However, the numerical result shows 
that they are well separated even at the lowest accessi- 
ble temperatures. This is a consequence of the different 
logarithmic corrections stated above. Quantitatively we 
can regard a multiplicative logarithmic term In" r as an 
effective, distance dependent correction of the scaling di- 
mension: 



1 ln(ln" r) 

2 Inr 



(18) 



The relevant length scale is given by the correlation 
length at the considered temperature, r w £.{T). Using 
this correction together with Eq. ( |l6| ) leads to splitting 
of the 2kp s-s and d-d CL at finite temperature and to 
an excellent agreement with the numerical results (see 
triangle downs in Fig ^ . 

Next, we regard the singlet pair CL's shown in Fig. ||. 
The CL leading over the entire temperature range is non- 




Eq. (12d) identifies this contribution as the 2kp part. 
The form of the algebraically decaying singlet pair corre- 
lation Ps (r, 0) at zero temperature is given by CFT 



P,{r, 0) = Cor-^" + Car-'" cos{2kpr) 



(19a) 



The non-oscillating term is due to the excitation 
(2, 1, ±1/2, =f1, 0, 0) leading to an exponent 



1 



edQ) 



(19b) 



and the 2kp part due to the excitation (2, 1, ±1/2, 0, 0, 0) 
leading to 



' eAQ) 2 • 

Thus, the coefficient 7 in Eq. (|l^ ) is now 

Vc 

7 ■ 



Jim ' V, 

for the non-oscillating part, whereas 



7 = 



'( 



;g(Q) 



11(0) 



(19c) 



(20a) 



(20b) 



for the 2kp part. The value for the non-oscillating part 
is again in good agreement with the numerical results 
whereas the numerical accessible temperature range is 
not sufficient to compare the 2kp part with CFT (see 
symbols in Fig. ^). 

To illustrate the dependence of the leading s-s and 
singlet-pair CL on electron density, there are plots for 
three different chemical potentials in Fig. |^. The s-s 



0.75 



0.5 



0.25 



— H=-0. 7 spin-spin 


1 ' 1^ ' 


— H=-0. 7 singlet pair 




--- (i=-1. 4 spin-spin 




-— n=-1. 4 singlet pair 




- - n=-1 .9 spin-spin 




- - n=-1 .9 singlet pair 








1 


1,1 



0.5 



1 

T/t 



1.5 



FIG. 6: Leading singlet pair correlation lengths for J/t = 2.0 
and ^ = —1.4. The square at T = denotes the CFT result 
for the non-oscillating part, the circle that for the 2kF part. 
The inset shows the corresponding wavevectors. 

oscillating, whereas the next-leading shows incommensu- 
rate oscillations. Again the zero temperature result from 



FIG. 7: The three upper curves show the spin correla- 
tion lengths, the other three the singlet pair-pair correlation 
lengths. 

CL is largest for the Heisenberg case {n — \) and gets 
suppressed, as expected, with an increasing concentra- 
tion of holes. On the other hand the singlet-pair CL is 
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nearly independent of the particle density. The s-s CL 
shows a non-analyticity at a temperature T^, where the 
k — TT part is crossed by the non-oscillating part. With 
increasing particle density this non-analyticity is shifted 
to lower temperatures and disappears for chemical poten- 
tials /i > —0.47. Therfore no singularity is visible in the 
Heisenberg limit of the t-J model. Tc as a function of the 
chemical potential shows an algebraic behavior. A very 
similar phenomena is described in Ref. |l^ where the s-s 
CL of the 1/2-XXZ chain in a magnetic field is investi- 
gated. In this case the magnetic field plays the same role 
as the chemical potential here. However, the singularity 
there is associated with a crossover from commensurate 
oscillations at T > Tc to incommensurate oscillations at 
T < Tc- Recently, temperature dependent CL's have also 
been studied analytically in a generalized Hubbard model 
related to the supersymmetric t-J modelJlj 

B. Static correlation functions 

Within the TMRG imaginary time correlation func- 
tions are easy to calculate.llj The fatal point is that the 
analytical continuation of the imaginary time results to 
real times is an ill-posed problem. We therefore want to 
restrict ourselves to the calculation of static correlation 
functions defined by 

G{r,z = 0)= / dTG{r,T) (21) 

where G{r, r) is the correlation function for distance r 
and imaginary time r. Note, that the static correla- 
tion function is calculated to higher accuracy than the 
equal-time correlation. Nevertheless the spatial CL's are 
identical in both cases except for the trivial correlations 
caused by particle-hole excitations which have vanishing 
matrixelements in the static case (see Eq. ^5|) . The static 
correlation function can be expressed through the largest 
eigenvalue and corresponding eigenvectors of the QTM: 

G{r, z = 0) = {^^iTMn^'TM^) (22) 

for distances r > I with the definition 

M 

Tm = Y,TM{A,.k) (23) 

where TM{A^-k) denotes the usual transfer matrix Tm 
where the considered operator A is added at imaginary 
time position t = e ■ k and 

G(r = 0,z = 0) = -^(vl/o^|fM|0 (24) 

with 

M 

Tm -^TM(Ao,A,.fc) . (25) 

fc=0 



These equations are valid within the alternative mapping 
to the classical model with alternating rowstj where a 
local transfer matrix covers only one site (see Sec. |). 
Note, that these formulas become more complicated if 
the usual decomposition of the Hamiltonian into even 
and odd parts is applied, because in this case even 
and odd distances as well as distance r = 1 have to 
be treated separately. To check the accuracy of the 
TMRG algorithm we have calculated the static density- 
density correlation function for free spinless fermions (see 
appendix |^ showing that the TMRG results are reli- 
able even at short distances r. In Fig. || (Fig. ^ the 




FIG. 8: Static longitudinal s-s correlation function at T = 2.0 
showing TT-oscillations. The dotted line denotes an envelope 
corresponding to the exponential decay. The solid line is a 
guide to the eye. 

static longitudinal s-s correlation function for J/t = 2.0, 
^ = -1.4 and T = 2.0 (T = 0.1) is shown. From 
Fig. ^ we expect 7r-oscillations of the correlation func- 
tion for T — 2.0 and 2fci?-oscillations for T = 0.1 be- 
cause the non-oscillating part does not show up in a 
static correlation function as explained before. By us- 




r 



FIG. 9: Static longitudinal s-s correlation function at T = 0.1 
showing incommensurate oscillations. The dotted line denotes 
an envelope corresponding to the exponential decay. The solid 
line is a guide to the eye. 

ing (SqS^) — Aexp(— r/^) cos(fcr + (5) as a fit function 
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we have estimated the CL's as well as the wavevectors 
for both temperatures directly from the correlation func- 
tions (see legends in Fig. ^ and Fig. ||). In both cases the 
results coincide within errors of the order 10"'' with that 
obtained by using directly the eigenvalues of the QTM. 
The perfect coincidence of the CL's estimated directly 
from the eigenvalues and that from a fit of the correlation 
function is not astonishing. In both cases the same QTM 
is used in the numerics where several eigenvalues are cal- 
culated directly using diagonalization routines in the first 
case whereas only the largest eigenvalue and the corre- 
sponding eigenvectors are needed to calculate the cor- 
relation function by some time-consuming matrix-vector 
multiplications in the second case. However, it provides 
a good consistency check of the numerics. 



IV. THE PHYSICALLY RELEVANT REGION 

As mentioned in the introduction the t-J model is a 
simplification of the Hubbard model obtained in the limit 
U ^ t. Thus, one is especially interested in the case 
where the spin exchange coupling J — 2t^ /\U\ is smaller 
than the hopping t. The TMRG is applicable for any 
value of J /t with almost the same accuracy and our algo- 
rithm remains stable even in the limiting cases J/t^ 1 or 
J/t 3> 1. Problems as reported by Ammon et alrB never 
occurred in our algorithm. To be specific, we never ob- 
served that the algorithm breaks down after performing 
only a few or more RG-steps although we have retained 
up to 240 states to calculate the correlation lengths ac- 
curately. In the following calculations we set J/t — 0.35, 
a value often used in the literature. 




FIG. 10: Specific heat for J/t — 0.35 and three different chem- 
ical potentials corresponding in the low temperature limit to 
a high (/I — 1.2), medium = —0.7) and low (/i = —1.4) 
electron density. Note again the two-peak structure due to 
spin-charge separation. The inset shows the temperature de- 
pendence of the densities n{fi, T) for the same chemical po- 
tentials. 

When regarding the specific heat shown in Fig. |l^ 
a two-peak structure is again obvious in the high and 



medium density case and there seems to be a linear 
regime consistent with TLL theory, although we have 
noticed in the supersymmetric case that numerical errors 
grow up for low temperatures. The structure of the spin 
and charge susceptibilities (see Fig. ^T]) looks also rather 
similar to the supersymmetric case. The arguments given 



X 1 -r 




FIG. 11: Susceptibilities and compressibilities for the three 
different chemical potentials and J/t — 0.35. 

for the high and low temperature regime of the suscepti- 
bilities are also applicable here. It is easy to understand 
that the absolute values of the spin susceptibility are now 
larger, because the weaker antiferromagnetic coupling J 
suppresses less efficiently a ferromagnetic alignment of 
spins. In principle, the qualitative discussion of c, Xs and 
Xc is very similar to that of the integrable case. There 
are, however, two modifications. First, in the charge sus- 
ceptibility for low densities we observe two maxima. We 
do not have any simple explanation of this fact. Note 
that a similar observation has been made for the Hub- 
bard chain. Second, due to the different energy scales of 
spinon and incoherent single particle motion at any den- 
sity there are two maxima in c(T) even at densities close 
to 1. 

We also calculated again several s-s and d-d CL's as 
shown in Fig. |l^. The CL's are now often difficult 
to distinguish and the crossing of CL's is "smoother" 
and shifted to higher temperatures when compared with 
the supersymmetric case. This clearly shows that the 
spin-exchange interaction is responsible for the various 
crossovers of CL's and in Appendix ^ we show that all 
such crossovers vanish for free spinless fcrmions. It is 
also striking that incommensurate oscillations are now 
present even at high temperatures. The leadin g s-s CL 
shows oscillations which are identified by Eq. (12d) as 
2kp whereas one of the d-d CL shows Akp oscillations. 
The other d-d CL from Fig. ( p^ ) is non-oscillating and 
leading in the low-temperature limit. These kind of os- 



cillations are expected from TLL theory (see Eq. (13a)) 



and we also note that the CL's seem to diverge again as 
^ ~ 1/T for low temperatures. This provides additional 
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- 0.3 



density-density (k-41(p) 

density-density (lc-0) 
spin-spin (Ic— 2ltp) — 2 

I I I I I L 



FIG. 12: S-s and d-d CL's for J/t = 0.35 and /i = -0.7. 
The inset shows the corresponding wavevectors. The circles 
(diamonds) denote the CFT result for ±2kF (±4^^). 



evidence that the t-J model at J/t = 0.35 belongs to 
the same universality class (TLL) as the supersymmetric 
model. 



V. PHASE SEPARATION 

At J/t large, the attractive Heisenberg interaction in 
Eq. (|l|) dominates the kinetic energy term. In a canoni- 
cal ensemble the model therefore phase separates into a 
high density and a low density region in order to optimize 
the Heisenberg energy. In our grand-canonical descrip- 
tion of the model, however, there remains a competition 
between the chemical potential term and the Heisenberg 
exchange energy. If we ignore the kinetic energy com- 
pletely - what is exact if total phase separation occurs - 
a simple picture evolves. Because the ground state en- 
ergy of a s = 1/2-Heisenberg chain per particle is given 
by — J In 2, we expect in the limit T — ^ an empty state 
if /i < — J In 2 and a state with n = 1 if fi > — J In 2. 
Therefore phase separation can only be present at one 
spec ial po i nt characteri zed by w — Jln2. In the Fig- 
ures 



13(a)| , |l3(b)| , |l3(c)| the density at a constant temper- 



ature is shown as a function of the c hemica l potential for 
three different parameters J/t. Fig. |l3(a) shows that at 
J/t — 3.1 the compressibility dn/d(x in the limit T ^ 
is not diverging, indicating that the phase se parate d re- 
gion is not yet reached. At J/t = 3.3 (Fig. 13(b) ) the 
density jumps from n = to n w 0.8 at /i w —J In 2 
in the limit T as expected for the phase separated 
region. The phase separation is here between the empty 
and an electron rich state. At J/t = 3.5 we find full 
phase separation, indicated by the ju mp of the density 
from to 1 for r — > (see Fig. 13(c) ). Our calculations 
confirm that the fully phase-separated state is destroyed 
by introducing holes into the Heisenberg chain island as 
stated by Ogata et alx 




/// 



(a) J/t = 3.1 



(b) J/t = 3.3 
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(c) J/t = 3.5 



FIG. 13: The figures show the density as a function of the 
chemical potential for fixed temperatures T in units of t. Note 
that the lines connect a finite number of data points. 



VI. LUTHER-EMERY PHASE 

As already mentioned in the introduction a phase with 
a spin gap is expected at least for low densities and val- 
ues of J/t close to phase separation. Because we expect 
the spin gap to be caused by a marginal operator leading 
to an exponentially small gap, the TMRG is not suited 
to determine the phase boundaries. However, in regions 
where the spin gap is larger than the lowest accessible 
temperature the TMRG can show the existence of such a 
phase without using any additional assumptions. A spin 
gap is directly visible in the spin susceptibility going to 
for T — > 0. Therefore, we have calculated the spin sus- 
ceptibility for two different values of J/t close to phase 
separation and have chosen chemical potentials so that 
the density is given by n « 0.2 for low temperatures (see 
Fig. lU). In both cases a very small spin gap appears. 
The quadratic dispersion of a gapped ID system leads 
to X '--^ exp(— A/r)/Vr for the low-temperature asymp- 
totics. Using this function for a fit of the numerical data. 
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J/t=2.9; \i=-2J3 
J/t=3.1;n=-2.2 




FIG. 14: Spin susceptibilities for two different parameter sets 
(J/t,^), both showing a spin gap of the order A ~ 0.05. The 
inset shows the temperature dependence of the corresponding 
densities. 



we find A = 0.05 ± 0.01 in both cases. 



0.4 



density-density (k=0) 
density-density (k=:7i) 
spin-spin (k=7r) 
spin-spin (k=0) 



1.5 
T/t 



FIG. 15: Correlation lengths at J/t = 2.9 and /i = -2.13. 
Note, that the oscillations of all shown correlation lengths are 
commensurate over the entire temperature range. 



Another proof of LE-properties of the t-J model is 
given by the calculation of s-s and d-d CL's. In Fig. |l^ 
these CL's are shown as usual as temperature times CL 
versus temperature. A spin gap is connected with fi- 
nite s-s CL's and we therefore expect that T • ^ 
for T ^ 0. On the other hand the d-d CL's £,n are not 
affected and should still diverge as ^„ ~ 1/T for low tem- 
peratures. This picture seems to be consistent with the 
numerical results. However, we are not able to present 
numerical data for lower temperatures, which could sup- 
port this scenario further. 



VII. T-J MODEL WITH ISING ANISOTROPY 

In this section we want to consider the t-J model with 
modification of Hamiltonian by replacing 



QX QX 



i+1 



(26) 



where A > 1. We expect that such an anisotropy en- 
hances superconducting correlations relative to d-d cor- 
relations as has been explicitly shown in an exactly solv- 
able anisotropic t-J model by Bariev et aLc3 However, 
their version includes some unphysical parity breaking 
terms (keeping PT-symmetry) making it interesting to 
investigate if the same is true for the model defined here. 
For the Heisenberg chain it is known that an Ising-like 
anisotropy promotes long range spin order and causes a 
spin gap. The situation is more complicated in the t-J 
model: The charge sector is unaffected and the charge 
excitations remain critical (i.e. d-d correlations decay al- 
gebraically in the ground state) . The expected long range 
spin order could be hidden when the density n 7^ 1, 
because the spin operators act on the physical lattice 
whereas the spins are coupled to the "electron lattice". 
The long range order would then be visible only in a 
string order parameter (SqS^) where all empty sites be- 
tween and r are omitted|^This has been emphasized 
also by Pruschke and Shibacll who have studied the limit 
J/t 0. In any case there has to be a spin gap also 
in the t-J model if A > 1. However, in the Heisen- 
berg chain the gap is given by AE oc exp {—J/ (A — 1)} 
leading to an exponential small gap and this gap is fur- 
ther reduced in the t-J model for densities n < 1 mak- 
ing it often undetectable in the thermodynamical data. 
In Fig. no a spin gap As — 0.054 ± 0.01 is visible in 




FIG. 16: Susceptibilities (main figure) and compressibilities 
(inset) for three different chemical potentials where J/t — 2.0 
and the anisotropy parameter is given by A = 1.5. 

the susceptibility data for /i = —1.5 corresponding to 
n « 0.98 in the low-temperature limit. This is in agree- 
ment with the gap A^^ = 0.043 for the Heisenberg chain 
with the same anisotropy. For ^ = —1.8 (tit^o ~ 0.54) 
and = —2.0 (riT-»o ~ 0.27) the spin gap is so small 
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that it is not visible in the accessible temperature range. 
The specific heat (Fig. |lj) looks qualitatively similar to 




^= 


1.5 




1.8 


^= 


2.0 



1.5 

T/t 



FIG. 17: Specific heats for the same chemical potentials and 
same parameters as in Fig. 

the isotropic case. Again two peaks (shoulders) are vis- 
ible corresponding to spinon and holon excitations. In 



long, spin-spin (k=jl) 
long, spin-spin {k=0) 
long, spin-spin (k=n) 
density-density (k=0) 
density-density (k-7t) 
singlet pair (k=0) 



2 

T/t 



FIG. 18; Leading d-d and longitudinal s-s CL's for /i = —1.8 
and A = 1.5. Additionally, the leading singlet pair CL is 
plotted. The shown CL's have commensurate oscillations. 

Fig. ^ several crossovers in the CL's are visible. In the 
low-temperature regime the leading d-d and longitudi- 
nal s-s CL's diverge as £, ^ 1/T showing that the charge 
excitations remain critical and that the long-range spin 
order is hidden. As expected the singlet pair CL dom- 
inates at low temperatures over the d-d CL supporting 
that an Ising-like anisotropy could be a mechauism lead- 
ing to superconductivity in higher dimensions.E3 



relation functions much easier. At the supersymmetric 
point of the one-dimensional t-J model an excellent agree- 
ment between the numerical results and Bethe ansatz was 
found. Additionally, the numerical results for correlation 
lengths and static correlation functions have been suc- 
cessfully tested for free spinless fermions. This shows 
that the TMRG is well suited even for fermionic sys- 
tems and yields accurate results over a wide tempera- 
ture range what was sometimes questioned before.c2l In 
particular, we have concentrated on the calculation of 
correlation lengths at the supersymmetric point to study 
the crossover from the non-universal high T lattice into 
the quantum critical regime (T <^ t). The non- universal 
regime is characterized by various crossovers between 
CL's with different wavevectors whereas the CL's are 
non-crossing and diverging as 1/T in the universal TLL 
regime. A good coincidence between predictions for the 
low-temperature asymptotics of CL's by CFT and the 
numerical results was shown, but it was important to 
take also the logarithmic corrections into account. For 
J/t = 0.35, a value often considered physically relevant, 
we observed properties that were rather similar to the 
supersymmetric case. In particular, the t-J model at 
this parameter point belongs also to the TLL universality 
class. In our grand-canonical description phase separa- 
tion has the meaning of phase coexistence at one chemical 
potential fi. This is due to the remaining competition be- 
tween the Heisenberg and the chemical potential term. If 
J/t is large enough so that total phase separation occurs 
this competition leads to an empty state for T = if 
/i < —J In 2 and to a state with n = 1 if /i > —J In 2. 
Our data also show that the fully phase separated state 
is destroyed by introducing holes into the Heisenberg is- 
land. By calculating directly the spin susceptibility for 
small densities and values of J/t near phase separation, 
we have proven the existence of a spin-gap phase without 
making additional assumptions. This was also supported 
by the calculation of s-s and d-d CL's, indicating that all 
spin CL's are finite at zero temperature, whereas the d-d 
CL's are unaffected and diverging still as 1/T. Finally we 
have studied the t-J model with an additional Ising-like 
anisotropy and have shown that singlet pair correlations 
are enhanced so that a tendency towards superconduc- 
tivity is expected in higher dimensions. Further on, we 
have settled that the expected long range spin order is 
hidden away from half-filling, because the spins are cou- 
pled to the "electron lattice" whereas the operators act 
on the physical lattice. 



VIII. CONCLUSIONS 
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APPENDIX A: FREE SPINLESS FERMIONS 

Here we calculate the density-density correlation func- 
tion for a system of free spinless fermions on a lattice. 
The Hamiltonian is given by 



H = —t cos(g)c^Cg 



(Al) 



and the density operator by 

^^(*) = -^E4+,/2Wcfc-9/2W- (A2) 

k 

By applying Wick's theorem the d-d correlation function 
separates into two-point functions 

t(i)c.(i)4(0)co(0)) = (^"^''^)) 
i E ^'"^ (4- f (*)^.-f (0)) (c..+ | (t)cl+, (0) 



k,q 



where = — icos(A:) and n(efe) is the Fermi function. 
By using the Matsubara formalism, inserting the known 
results for the electrons Green's function and then going 
back to real time, this is easily transformed into 



4{t)cr{t)cl{0)co{0)) = 

\ piqr 

q^k 



(A4) 



l)(e 



1) 



Regarding only the connected part and transforming to 
Matsubara frequencies we get in the free fermion case 



G{r,z = 0) = 



1 



dq dk — cos(gr) 



(A5) 

for the static density-densitycorrelation function as de- 
fined in Eq. (|||). In Fig. ll| we show results for two 
characteristic temperatures obtained by TMRG using 
Eqs. (IH) and in comparison t o th e exact results 
from a numerical evaluation of Eq. (A5). The absolute 
error is largest for short distances and decreases exponen- 
tially with distance. However, even at T/t — 0.1 the error 
always remains smaller than 10"^ showing that TMRG 
yields useful results not only for the asymptotic behavior 
of static correlation functions but also at short distances. 
Returning to Eq. (A4) and setting t = we find 



1 r e"?'' 

2tT / ^ '^e~Pcos(q) 



1 

2^ 



-ikr 



dk- 



.J3 cos(fc) ^ 



(A6) 



Now we can in principle evaluate these integrals by clos- 
ing the integration path in the upper (lower) half plane 




BO 



I 



|o]0 4 8 12 16 20 24 28 32 



20 



32 



FIG. 19; The upper (lower) graph shows the static density- 
density correlation function for fre e fermions at T/t = 1.0 
{T/t = 0.1) calculated by Eq. (^ (circles) and by TMRG 



(squares) with A'^ = 120 states retained. The correlation 
function shows in both cases commensurate oscillations with 
k — TV. The insets show the absolute errors (triangles) and 
the relative errors (diamonds) of the TMRG results. The lines 
are guides to the eye. 



and using Cauchy's formula. The poles are at the points 



TT / TT 

g, fc = ± — ± i arcsinh (2?! + l)-^ 



and thus, the dominant contribution is 

(^clcrclco) - Me-'-/«(l - cos(7rr)) 

where the matrixelement is given by 

2 



M 



(3^ 



and the CL by 



2 arcsinh(7r//3) 



(A7) 



(A8a) 



(A8b) 



(A8c) 



In Fig. ^ the leading CL received by Eq. (ASc) is shown 
together with TMRG resuhs calculated by Eq. (|). A 
good agreement is obtained with errors of the numerical 
data remaining smaller than 10"^ for temperatures down 
to T/t ~ 0.1. 
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